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Abstract 

We compute the two-point correlator between left-handed flavour charges, and the three-point 
correlator between two left-handed charges and one strangeness violating A/ = 3/2 weak 
operator, at next-to-leading order in finite volume SU(3) l x SU(3)/j chiral perturbation theory, 
in the so-called e-regime. Matching these results with the corresponding lattice measurements 
would in principle allow to extract the pion decay constant F, and the effective chiral theory 
parameter 52?) which determines the A/ = 3/2 amplitude of the weak decays K — > vrvr as 
well as the kaon mixing parameter Bk in the chiral limit. We repeat the calculations in the 
replica formulation of quenched chiral perturbation theory, finding only mild modifications. 
In particular, a properly chosen ratio of the three-point and two-point functions is shown to 
be identical in the full and quenched theories at this order. 
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1. Introduction 



The study of the weak matrix elements involved in kaon physies is a long-standing topic of 
lattice QCD [1, 2]. It is however a very difficult problem for all lattice formulations which 
break the chiral symmetry explicitly, such as Wilson fermions. We expect that significant 
progress can be achieved with the new Ginsparg- Wilson formulations of lattice fermions [3]- 
[10], which possess an exact chiral symmetry in the limit of vanishing quark masses, resulting 
in an enormous simplification in weak operator mixing and renormalization [8], [11]-[13]. 

One advantage of approaching the regime of vanishing quark masses is obviously that the 
uncertainty induced by chiral extrapolations is avoided. On the other hand, since the chiral 
limit implies that light mesons become massless, it necessarily brings with it large finite 
volume effects. This apparent difficulty can, however, be turned into a useful tool: if the 
finite volume effects can be resolved analytically in terms of the infinite volume properties of 
the theory, then the infinite volume properties can be extracted by monitoring the volume 
dependence. The first proposal to apply finite-size scaling techniques to the weak K mr 
amplitudes was presented in [14]. 

A practical realisation for the finite-volume philosophy mentioned is offered by Chiral 
Perturbation Theory (xPT). As the quark masses get smaller, the chiral expansion becomes 
more and more accurate in describing the dynamics of the low momentum modes of QCD 
(below a few hundred MeV). The chiral expansion in this regime is slightly more complicated 
than in infinite volume, because it requires the resummation of pion zero mode contributions. 
Gasser and Leutwyler [15] have presented a systematic procedure for doing this, the so-called 
e-expansion (see also [16]). Several observables, such as the quark condensate and the scalar 
and vector two-point functions, have already been computed at next-to-leading order in 
the e-expansion [15]-[18]. These quantities depend on the (infinite volume) chiral theory 
parameters, the volume, and the quark masses, in a way that a comparison with lattice data 
for different volumes and quark masses, allows in principle the extraction of the corresponding 
infinite volume low-energy couplings of xPT. 

Obviously the chiral model can be extended to include the \AS\ = 1 weak Hamiltonian [19]. 
This introduces a new set of low-energy constants, from which the physical amplitudes in 
kaon decays can be extracted, by working up to some desired order in the chiral expansion. 
The determination of these low-energy constants by matching the matrix elements of weak 
operators computed in lattice QCD to the same observables computed in SU(3)LxSU(3)fl 
xPT was proposed a long time ago [2]. Calculations along these lines (for a recent review, 
see [20]) have shown that if the matching is performed at relative large quark masses, there 
are large uncertainties induced by the chiral extrapolations [21], due to the fact that next-to- 
leading order corrections involve a large number of new unknown couplings [22]. Although 
strategies have been proposed [23] to measure the relevant new couplings, together with the 
leading order ones, by matching the matrix elements at several kinematical conditions, this 
is clearly a very challenging procedure, particularly for the AI = 1/2 kaon decays [23, 24]. 
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The approach that we consider here is instead to perform the matching in a finite volume 
but close to the chiral limit, in the e-regime. The predictions of xPT for the weak matrix 
elements in terms of the low-energy couplings are then not the same as in infinite volume. We 
compute the correlators of two left-handed flavour currents, as well as the matrix elements 
of the IAS"! = 1 (or \AS\ = 2) weak Hamiltonian with two such currents, at next-to- leading 
order in the e-expansion. Some motivations for this approach have been discussed in [13]. 
The simplification brought in by approaching the chiral limit is manifest in the fact that at 
next-to-leading order none of the unknown couplings present in the usual p-expansion [22] 
contributes. We expect therefore that the determination of the leading order couplings should 
in principle be more straightforward. Only the 27-plet low-energy coupling contributing to 
AI = 3/2 kaon decays is considered here, due to a number of subtleties with the octet 
operator [23, 24], which will be considered elsewhere. 

The structure of the paper is as follows. In Sec. 2, we discuss the chiral model representation 
of the weak Hamiltonian at low energies in the SU(3)LxSU(3)ii symmetric case. In Sec. 3, 
we review the e-expansion of Gasser and Leutwyler for this model, and in Sec. 4 present 
the results of our ncxt-to-leading order calculations. Finally, in Sec. 5, we present the same 
results in the quenched theory. We conclude in Sec. 6. 



2. Weeik operators in the SU(3) chiral theory 



Ignoring weak interactions, the QCD chiral Lagrangian possesses an SU(3)lxSU(3)j? sym- 
metry, broken "softly" by the mass terms. The Euclidean Lagrangian can to leading order 
in a momentum expansion be written as 



-Tr 



-Tr 



(2.1) 



Here U G SU(3), 9 is the vacuum angle, Nf = 3, M is the quark mass matrix and, to leading 
order in the chiral expansion, F, S, equal the pseudoscalar decay constant and the chiral 
condensate, respectively. We shall for convenience take M to be real and diagonal. 

Weak interactions break explicitly the SU(3)lxSU(3)r symmetry of Eq. (2.1). In the 
fundamental theory, the strangeness violating interactions responsible for kaon decays can be 
accurately described through an operator product expansion in the inverse W-boson mass. 
In the CP conserving case of two generations, the effective weak Hamiltonian to leading order 
in the QCD coupling constant is then (see, e.g., [25, 26]) 



+ R.C. 



(2.2) 



where Gp is the Fermi constant, Vij arc elements of the CKM-matrix, and P± = (1 it 75)/2. 
Taking into account QCD radiative corrections, the coefficients get modified, but Hw can still 
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be written as [27, 12] 

= 2V2GFVudVd E KiPwYsuud - PwYsccd) + hm[Om]sd] + H.c. , (2.3) 



where h^,hm are rcgularisation dependent dimensionless Wilson coefficients, and we have 
introduced the notation 

[Ou,]rsuv = {A7^.P-'^u){1psl^^P-A) , (2-5) 
[Omld = {ml-ml){[^pM]sP-^d + i^sP+[MH]d}. (2.6) 

Here r, s, u, v arc generic flavour indices, while u, d, s, c denote the physical flavours. The 
Wilson coefficients have been computed also for Ginsparg- Wilson "overlap" fermions [12], 
apart from hm, which remains undetermined. According to Eq. (2.2), the leading order 
values are = I, hm = 0. 

In order to match the Hamiltonian of Eq. (2.3) to the one in the SU(3) chiral theory, the 
first step is to decompose it into irreducible representations of the SU(3)LxSU(3)ij flavour 
group, present at low energies. For completeness, we review the general formulae for the 
decomposition in Appendix A. The weak operators are singlets under SU(3)ij, and denoting 
projected operators transforming under representations of SU(3)i: with dimensions 27, 8 by 
[Ow]tsuvj [Pw]ru' respectively, the weak Hamiltonian can be rewritten as 

= 2V2GFV^,V:,!^ht,[dXuud + lh+[RX,-h-[R^]:, 

- \{ht, + h-)[Oyj]sccd-^{ht-h-)[Oyj]scd^ + hrn[Orn]sd^ +ii-C. , (2.7) 

where 

1 1 

Pw]tuud = 2!^^"']*""'^ t^"']*"'^" ~ 5 (iQwjskdk + [Owjskkd^^ , (2-8) 

k=u,d,s 

skdk ± [Ow]skkd^ ■ (2.9) 

k=u,d,s 

The first operator in Eq. (2.7) transforms under tlie 27-plet of the SU(3)l subgroup: it is 
symmetric under the interchange of quark or antiquark indices, and traceless. The remaining 
ones, transforming as 3* (g) 3 and being traceless, belong to irreducible representations of 
dimension 8. 

The next step is to find the chiral analogue for this weak Hamiltonian, as well as for the 
left-handed fiavour currents, which we will use as external probes. In a convention for fermion 
fields where the Euclidean Lagrangian reads (7^^ = 7^, {7/^,7;^} = 25 ^y) 

Le = i>{l^^D^ + MP_ + M^P+)i^, (2.10) 
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we may define a left-handed current as 

(•^^)qCD ^ ii'rTrl^^P-i^u, (2-11) 



where the are Hermitean generators of the flavour SU(3)^. As usual, it is convenient to 
introduce an external left-handed flavour gauge field source, A^^, such that 



M = (2-12) 

Observables including (J^)qcd can then be addressed within the chiral theory by coupling 
also the pion field covariantly to A^, and taking functional derivatives with respect to it. 
More concretely, the partial derivatives of Eq. (2.1) are promoted to covariant ones, 

d^U ^ B^U = [d^ + iA^^T'^p, (2.13) 

and the left-handed current is defined as 

M Af=0 

up to higher order corrections. 

We can now find the chiral analogues for the building blocks in Eqs. (2.4)-(2.6). By a 
comparison of Eqs. (2.5), (2.11), the analogue of the operator in Eq. (2.5), denoted in the 
chiral case by [O^lrswu) can be written as 

[0^]rsuv = \F^id,/JU^) id^UU^) . (2.15) 

The operator of Eq. (2.15) is the only one with the same symmetry properties as its counter- 
part in the fundamental theory, at the leading order in the chiral expansion. Correspondingly, 
we can also find the chiral counterpart for the operator O^, by employing scalar and pseu- 
doscalar external sources S"", P", defined with the substitution M M + S"'T"' — iP"'T°- , and 
taking derivatives with respect to the sources: 

PmU = -{ml - ml)^{UMe''l''f + M^U^ e-'^'^^nds • (2-16) 

Given these building blocks, we can determine all the operators (at the leading order in the 
chiral expansion) transforming under the 27-plct and octet of SU(3)l, which allows then to 
translate the weak Hamiltonian of Eq. (2.7) to the chiral theory. Because some contractions 
are traceless, cf. Eqs. (A.8), (A. 9), there are only three such operators [2]. For the physical 

^In fact, the only property of T" we need to assume is their tracelessness. 
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choice of indices, we write these as 

O27 = [<^w]tuud = ^{P^]sudu + ^P^flsuud)^ (2-17) 

k=u,d,s 

O's = ^F^^{UMe'^/^f +M^U^e-'^/^f)ds , (2.19) 

where we have made use of Eqs. (A. 8), (A. 9) to simphfy the chiral versions of Eqs. (2.8), (2.9). 
Note that in the definition of Og here we have left out the explicit mass combination (m^—m^), 
which can then appear in the coefficient of this operator; the coefHcient can, however, also 
receive other contributions, due to mixings with operators of the same symmetries. 

We can now write down the analogue of in Eq. (2.7) in the chiral theory. We denote it 
by Hw To again define dimensionless coefficients, we write Hw in the form 

Hy, = 2V2GFVudV:, I ^527^27 + 258^8 + 2g'sO's | + H.c. , (2.20) 

where 527,58 and g'g are the low-energy constants we are interested in [28]. 

Now, it is easy to sec that the amplitude for A/ = 3/2 decays, such as — ^ tt^tt'^, is 
directly proportional to 527, while the much faster AI =1/2 decays of Kg get a comparable 
contribution both from and (727- (The parameter g^, on the other hand, does not contribute 
to physical kaon decays [27, 2, 29].) More quantitatively, a leading order analysis in infinite 
volume [28], supplemented by phenomenologically determined large phase shifts [30] in the 
amplitudes, suggests the well-known values 

Iff27| ~ 0.29, (2.21) 
|<78| « 5.1. (2.22) 

It has been argued that 1-loop corrections in the chiral perturbation theory are large [22, 
31, 32], and one can therefore get agreement with experimental data on partial decay widths 
even with somewhat less differing values of gs and 527, but a hierarchy still remains. 

In the limit iVp — 00, on the other hand, one obtains [28] the "tree-level" values deducible 
from the naive conversion of Eq. (2.7), with as in Eq. (2.2), to the corresponding chiral 
operators of Eqs. (2.17), (2.18): 

527 = 58 = ^ ■ (2.23) 

Clearly a large non-perturbative enhancement of gs with respect to the tree-level value, and 
some reduction of 527, is needed to fit the experiment. The final goal is to improve on the naive 
estimates in Eq. (2.23), by determining 527 and gs non-perturbatively in the SU(3)LxSU(3)i? 
symmetric theory. As mentioned in the introduction, we will in this paper discuss only 527, 
due to various subtleties in the determination of gs (particularly in the quenched case). 
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Let us finally recall that another physical observable determined by 527 is the Bk, charac- 
terising the mixing of and hence determining the mass difference of Ks,Kl [33]. It 
is defined by^ 

( I hAS=2 [O^lsdd I ) = lirriKFKfBK, (2.24) 

where ^a5=2 is the Wilson coefficient (see, e.g., [34] and references therein) related to the 
operator [Ow]ssdd, normalised to unity at tree-level. Since [Ou,]ssdd is symmetric and traceless, 
it belongs to the 27-plet. Therefore, if /iAS=2 is replaced by in Eq. (2.24), the matrix 
element is in the chiral limit (where uik = m-;^ = 0, Fk = F^^ = F) proportional to g2T- 

4 A hAs=2 5 

The tree- level value is then Bk = 3/4, but going to next-to-leading order in the large- A?c 
approach one finds a suppression down to Bk — 0.3... 0.4 [35]. This suppression factor is very 
close to what would be needed for (727 to go from Eq. (2.23) to Eq. (2.21). We may thus 
consider it a further motivation for the lattice study to corroborate this prediction of the 
large-ATp approach. Note that near the physical point (m^ > 0) a considerably larger value 
is found (for recent reviews, see [36]). 



3. Chiral perturbation theory in a finite volume 

In order to determine 527 , we will consider the lattice measurement of left-current two-point 
correlation functions in the fundamental theory at low enough momenta [13]. In this regime 
we assume that the effective theory gives a good description, and thus require, to first order 
in the weak Hamiltonian, 

where is set to zero after the differentiations. The measurement of the left-hand-side in 
lattice QCD allows to tune the effective couplings in the weak Hamiltonian appearing on 
the right-hand-side. In general, it is convenient even to remove the integral J^{---) appearing 
in Eq. (3.1), since matching can also be achieved before this averaging, as long as the currents 
brought down by the functional derivatives arc far enough from each other, and H^{z). In 
this way we avoid complications with "contact terms", arising from operators overlapping 
at the same spacetime location. It is also convenient to consider space-averaged charges 
positioned at different times, xo,yo, rather than local currents. 

We thus discuss the product of two left-handed charges separated from the weak Hamilto- 
nian sitting at the origin, z = 0. To keep the discussion as general as possible, we consider 

^The parameter Bk is defined identically to Bk but without the Wilson coefficient /ias=2, whereby it is 
scheme and scale dependent. 
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matrix elements of the "unprojected" operator [C>«,]rsuj; in Eq. (2.15); a projection to the 
actual 27-plet [Ojul^^^ is then carried out by the operations in Appendix A. Thus, writing 
the expressions in a form where their QCD analogues are obvious^, we will be concerned with 

C^^xo) = [{X{x)Jo'm, (3.2) 

[Cw]fsuvi^O,yo) =11 {JS{x)[0^]rsuvm^{y)), (3-3) 
Jx Jy 

where = J d^x. The computations are carried out with a finite phase 9 as in Eq. (2.1), 
allowing to make predictions for the case of a fixed topology, as well (see Sec. 4.3). 

Due to the vicinity of the chiral limit, we take the volume to be a finite periodic box, of 
size V = L0L1L2L3. Momenta are then quantised, 

Pij. = T'nij,, e Z. (3.4) 

Since we want to be close to the chiral limit in the finite volume, the computation is organised 
according to the rules of the e-expansion [15]. In the e-expansion one writes 

U = exp{i^)Uo, (3.5) 

where ^ has non-zero momentum modes only, while Uq is a constant SU(3) matrix collecting 
the zero modes. The integration over Uo has to be carried out exactly when m'EV <0{l), 
where m is a quark mass, while the integration over the non-zero modes can be carried out 
perturbatively as long as FL S> 1. The power counting rules for the e-expansion are 

F^Oil), 9^~0(e), L^~0(l/e), ^ ~ 0(e), m ~ ^(e^). (3.6) 

Note that the quark mass counts as four powers of the momenta, rather than two as in the 
standard chiral expansion in infinite volume, where m ~ 5^. The perturbative integrals 

for the non-zero momentum modes arc computed with the measure 

We will compute at next-to-leading order in the e-cxpansion, including corrections of rela- 
tive order 0{e^). It turns out that the physical pion decay constant and mass, F-j^jM-,^, differ 
from their leading order values, F, 2mS/F^ (for M = diag(m, m, m)), only by terms of rela- 
tive order O(e^) [18], an effect which may thus be ignored. This is a result of the fact that no 
higher order operators in the action (i.e., none of the Lj's of Gasser and Leutwyler) contribute 
at 0{e^). Similarly, there are no higher order operators of C(e^) in the chiral representation 

*At the present order, these forms differ from those obtained by taking functional derivatives with respect to 
a flavoured gauge fleld, as in Eq. (3.1), only regarding unimportant "contact terms" ~ 5{xo),5{yo), 5{xo — yo). 
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T/L 






32/32 


0.14046 


0.07023 = /3i/2 


32/28 


0.13872 


0.07826 


32/24 


0.13215 


0.08186 


32/20 


0.11689 


0.08307 


32/16 


0.08360 


0.08331 



Table 1: Some numerical values for fcoo, defined in Eqs. (3.10), (3.11), for geometries of 
the type Lq = T, Li = L2 = = L. 



of [Otujr-sui; in Eq. (2.15). This is also in contrast with the usual chiral expansion in infinite 
volume, where a plethora of new operators contribute at next-to-leading order [22]. The fact 
that the contamination from higher order operators is small, is simply a consequence of being 
closer to the chiral limit. 

The propagator for the non-zero momentum modes ^ follows by expanding the parametri- 
sation in Eq. (3.5) in e, and inserting into Eq. (2.1): 

Cur{x) ^vs{y)) = 7; UusSvrG{x - y) - 6ur6vsE{x - y)] , (3.8) 



2 

where 

G{x)= (3.9) 
Jp' p 

In the chiral case, E{x) = G{x)/Nf, but we keep everywhere E{x) completely general. The 
reason is that then the form of Eq. (3.8) is general enough to contain also the propagator 
of the replica formulation of quenched chiral perturbation theory [37, 38], and thus we can 
already include the main ingredients needed in Sec. 5. 

Where encountered, ultraviolet divergences are treated with dimensional regularization. 
For later reference let us define [39, 17], in particular, the two integrals appearing in the 
computation (see also [40]): 



1 _ Pi 



p2 - Y1/2-- 



(3.10) 



Here Pi , koQ arc finite dimcnsionlcss mimcrical coefficients depending on the geometry of the 
box. Some values for them are listed in Table 1. 
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□ — H — □ 



0. 



Figure 1: The O(e^) graphs computed in Sec. 4-1- An open box denotes a current (Eq. (2.14)), 
a cross a "measure term" (cf. ref. [15]), and a filled circle a mass insertion. 

4. Results 



4.1. Charge — charge expectation value 

We now proceed to apply the rules of the e-expansion to the correlator ^^^(xo), defined 
in Eq. (3.2). The result can in fact also be inferred from [17], by summing together the 
expressions for the axial and vector flavour currents. 

The graphs contributing to C"*(xo) are shown in Fig. 1. Apart from the graph including the 
mass insertion, current conservation guarantees that the result is independent of xq. Indeed, 
we obtain 



(-TrTT^ 



Nf 
H 



2Lo 

Pi 



V 



+ 



p|(ReTY[MC/oe"^/^/: 



e,Uo \L 



where /3i,A;oo stre from Eqs. (3.10), (3.11); and [39] 



h(r) = 5 



' 2/ 



1 

12 



Finally, (...)e,C7o denotes an average over the zero-mode Goldstone manifold, 

/^^(...) exp(ySReTr [MUoe''/^f]) 



{■■■)e,Uo 



exp(l/S Re Tr [MUoe^'^/^f]) 



(4.1) 



(4.2) 



(4.3) 



where Jjj^ is an integration over SU(A^/) according to the Haar measure. 

As a simple explicit example, let us assume a box with Lq = Li = L2 = L3, a mass matrix 
M = diag(m, m, m), and a phase ^ = 0. Taking furthermore into account that at the present 



order 0{e'^), = F, = 2mT,/F^, the fuh result may be expressed as 



C'^^(xo) = (-Trr'^r'')^|l + -J- 
where [15, 17, 18] 



F2r2 



-di+uT.e=Q{u/2)h_ 



(4.4) 



MlF^^Ll 



(4.5) 
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T = 32a o"' = 2GeV 



L=24a T = 32a a'' = 2GeV 



1,4 


- 


1,4 


- 


1.2 


-^L=24a 


1.2 




ti. 




Ll. 




\ 




\ 

f< 1 


----'^ -y^---- 


b 




b 




' 0,8 




' 0,8 








0,6 


, , , 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 , , , , 1 


0,6 





5 10 15 20 25 30 5 10 15 20 25 

t/a t/a 



Figure 2: Le/i; T/ie expression inside the curly brackets in Eq. (4-1), for Lq = T , V = TL^, 
= 2 GeV, and m = MeV (dashed), m = 5 MeV (dotted). We have assumed 6 = 0, 
F = 93 MeV, T, = (250 MeV)^. The solid line is the tree-level result. Right: The same 
observable, in an ensemble with a fixed topological charge v = 0,1 (cf. Sec. 4-3), for L = 24a. 

^ ' ^ Nf\ /e,Uo Nfdu Juo ' v / 

W,/2) . {«• "<<;. (4.7, 
I 1, ll > 1 

Note that in the notation of [18], uT,0=o{u/2) = h{u) / {AN f) . 

For Nf = 3, a numerical determination of Ii{u) has been given in [18]. Using this result, we 
show in Fig. 2 examples for two asymmetric lattices, 32 x 24"^, 32 x 20^, and a lattice spacing 
= 2 GeV. The function in Eq. (4.1) has been normalised to its (constant) tree-level value. 
The next-to-leading order correction is observed to become dangerously large for L/a< 20 
(i.e., L<2fm). 

4.2. Charge — charge — weak operator Ow 

We then move to [CTO]rs„„(a;0) 2/o)) defined in Eq. (3.3). To this end we evaluate the graphs 
in Fig. 3. Let us mention that we are ignoring disconnected diagrams, since they only lead 
to trace parts, ~ S'^^Sur^vs^ ^''^^^usSvr, which would vanish in any case after the projection to 
[Owj^gy^jj. The non-trivial flavour structures arising from Fig. 3 always appear in one of the 
two combinations, 

r A (l)ia6 _ rpa rph , rpa rpb — rp{arpb} (A o\ 

L Wsuv — ur vs ' vs ur — -^ur vsi V 
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□ — H- 



-o □ n m o □ □ ^ 




Figure 3: The O(e^) graphs computed in Sec. 4-2- An open circle denotes the weak operator 
in Eq. (2.15), otherwise the notation is as in Fig. 1. 



ab 

rsuv 



(4.9) 



After lengthy but straightforward algebra, we obtain^ 



[Cw]rLvi^O,yo) 



All 



Pi L^koo 



+A(i)|^/Re'IV[MC/oe^^/^/]^ 



V 



where we have for clarity omitted the indices from [A^^^]"^^^, [A^^^]^^^^. Apart from the 
graph including the mass insertion, current conservation guarantees again that the result is 
independent of xq, yo- 

Once the flavour structure is projected onto the 27-plet according to Appendix A, we get 
our final result, 



[C27]rL(^0,yo) 



mx)[dXsuvmoiy)) ■ 



X Jy 



It is directly obtained from Eq. (4.10), by simply making the replacements 



r A (i)ia6 



/\ab 7 — 12 

'-^rsuv ; ^ -"^5 ^ 



(4.11) 



(4.12) 



where 



A 



ab 



\ [t^sT^I + T^T^]) + ^ {5us6vr + 6ur5us) Tr T'^T' 



(4.13) 



As a simple explicit example, we again assume a box with Lq = Li = L2 = L3, a, mass 
matrix M = diag(m, m,m), a phase 6 = 0, Nf = 3, and choose generators such that 

T^d^uSjs, T^j = 5idSju. (4.14) 

^The full expression before the volume average /^(••■) can also be obtained from the authors on request. 
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T=32a o" = 2GeV 



L=24a T = 32a a'' = 2GeV 



..••'L=24a 



O 

•-0,8 



0,6 - 



,.-•••1=200 



10 15 20 25 30 

t/a 




1/ = 



Figure 4: Le/t; The result of Eq. (4-10) for the index choice in Eqs. (4-14), (4-15), normalised 
to the tree-level value (solid line), for the same parameters as in Fig. 2. The upper set is for 
L = 24a, the lower for L = 20a. The dashed and dotted lines correspond to m = MeV and 
m = h MeV, respectively. Right: The same observable, in an ensemble with a fixed topological 
charge i/ = 0, 1 (cf. Sec. 4-3), for L = 24a. 



We also choose the physical indices for [Oiu]^„„, according to Eq. (2.17). Then A"^^^^ = 2/5, 
and the full result is 



p27]w(^o,yo) = 



2 

-—^-i 1- 
4Lg5^' 



1 



(4.15) 



where the notation is as in Eq. (4.4). 

In Fig. 4 we show the predictions of Eq. (4.10) for this index choice, normalised to the tree- 
level value, for two asymmetric volumes. Since there is only one 27-plet operator in Eq. (2.7), 
a measurement of 527 can then be obtained through the matching of the chiral three-point 
function with the corresponding lattice QCD measurement: 



927[C27]sLdixo,yo) 



I I (Jo»(x)[Q 



^ i suud 



(4.16) 



4.3. Fixed topology 

The results above were obtained in a fixed ^-vacuum. We can also perform a Fourier transform 
in 6 to obtain averages in sectors of "fixed topology" v [41]. This is interesting because in 
the quenched theory we expect to find poles in m in fermion propagators, which become 
dominant in the e-regime when mTiV <C 1. In the ^-vacuum there are no such poles in 
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the full theory, because topological configurations are strongly suppressed by the fermion 
determinant. However, when considering averages in sectors of non-zero fixed topology, the 
same poles are expected to appear in the quenched and the full theories. It is quite remarkable 
that these poles appear also in the corresponding effective chiral theories! Even though their 
presence does not affect the counting rules of the e-expansion, because mT,V is formally 
counted as a quantity of they obviously modify the chiral limit. The question we want 

to address here is whether there are such poles in the observables of Eqs. (3.2), (3.3). 

An observable in the sector of topological charge u, f^^ can be obtained from the observable 
in a ^-vacuum, fg, by 

/. = ^y^ d^e-^Ve. (4.17) 
In particular, assuming again M = diag(m, m, m) and defining 

Ze{u/2)= i e("/2)R'^^[f^o^^PW^/)], Z^{u/2) = ^ f'^'' dee-'^^Ze{u/2), (4.18) 
JUo 27r Jo 

the combination appearing in Eqs. (4.1), (4.10) gets replaced as 

S,(V2) = ^|^lnZ,(u/2) _l|.inZ.(n/2)^S4V2)- (4.19) 

The function is known [42, 41] to be 

Z^{u/2) = det[h+j-i{u/2)l (4.20) 

where the determinant is taken over an Nf x matrix, whose matrix element (i,j) is the 
modified Bessel function 

Thus, at fixed topology the results corresponding to Eqs. (4.4), (4.15) are obtained by the 
substitution Sg(n/2) — S,/(n/2). For small and large u we have (independent of iV/), 

W2)»p|''y"' (4.21) 

1, u > 1 

As expected the low mass behaviour {u <C 1) is drastically modified with respect to that 
in Eq. (4.7). This implies that even though the correlators remain finite for m ^ (i.e., 
there are no poles, because Y,y{u/2) is multiplied by u), their time dependence does not 
vanish. This is illustrated for C"''(xo) in Fig. 2 and for '^iiYjiud^^^^V'^) ™ ^^S- 4. 

4.4. Normalised correlators 

The predictions of the previous sections depend at leading order on the chiral theory parame- 
ter F, show rather bad convergence at small volumes, L<2 fm, and have (for a non-vanishing 
mass, as well as for a non-zero topological charge) a non-trivial dependence on xo,yo- All 
these dependences are a nuisance for the determination of g27- Fortunately, there seems to be 
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Figure 5: The result of Eq. (^.22), for the same parameters as in Fig. 2. At this order, the 
outcome is independent of XQ,yo. 

a large cancellation if we normalise the three-point function [^27]^^^^ by two charge — charge 
correlators, at least at the present order in the e-expansion. 

More precisely, let us again choose the indices in Eqs. (4.14), (4.15), and denote by C""^ (xo) 
a charge - charge correlator obtained by using the generators T", {T"")^ in the currents. 
Expanding the denominators, we then obtain that to relative order C(e^), 



rr 1 (r Mg27]gLd(^0,yo) , 2 



yV2 V 



The same result holds at fixed topology, as defined in Sec. 4.3. 

Thus, the time, quark mass, and topology dependences cancel completely in the ratio 
of Eq. (4.22)! The next-to-leading order correction is also numerically smaller than the corre- 
sponding corrections in the numerator and denominator separately. This result is illustrated 
in Fig. 5 as a function of the spatial volume, for a fixed time-like extent T = 3.2 fm. 

To summarise, the optimal method for determining 5(27 would appear to be from the equal- 
ity 

_ 3^+ [C27]norm(xo,l/o) 
527 - ■^l^wJ^T-i 7 T > (4.23) 

5 [C27jnorm(a;o,yoj 

where [C27]norm is the QCD-corrcspondcnt for the expectation value in Eq. (4.22). The 
independence of the outcome on the volume, quark masses, xo,yo, and topological charge, 
serves as a test of whether the regime of applicability of Eq. (4.23) has been reached. 
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5. The quenched case 



5.1. The basic setup 

Due to the numerical cost of dynamical Ginsparg- Wilson fermions, practical lattice simula- 
tions will, for a while still, have to resort to the quenched approximation. It is therefore 
of interest to study how the results of the previous sections are expected to be affected by 
quenching. The tool for this is quenched chiral perturbation theory, applied to the e-regime. 
Previous results in this setup exist for the quark condensate [43], and the scalar and pseu- 
doscalar [38] as well as flavoured vector and axial-vector [44] two-point functions. 

There are two approaches to quenched QCD, believed to be equivalent: the so-called 
supersymmetric (SUSY) formulation [45, 46] and the so-called replica method [37, 38]. In 
the former, bosonic "ghost" quarks are introduced in order to cancel the effects of the physical 
quarks; in the latter, the computation is carried out by keeping separately track of the Ny 
"valence" quarks appearing in the external sources, and the Nf dynamical quarks, and the 
quenched limit is obtained by taking Nf ^ for a fixed Ny ^ 0. 

The two methods are formally equivalent at the quark level, however their low-energy 
effective theories appear to be quite different. Assuming that the naive chiral symmetries of 
these models, \J{Nf)£^x\J{Nf)ii in the replica case and the graded \J{Ny\Ny)i^x\J{Ny\Ny)ji 
in the SUSY formulation, are broken spontaneously by the formation of a quark condensate 
to the corresponding vector subgroups, the low-energy degrees of freedom are the resulting 
Goldstone bosons, whose dynamics can be described by chiral Lagrangians, at energies below 
the typical confinement scale. The field variables of the chiral Lagrangians are matrices 
parametrising the Goldstone manifolds. 

There is one important difference with respect to full QCD, though: the field associated 
with the singlet axial rotation (~ r]') does not decouple from the low-energy dynamics. This 
is true both in the SUSY method [45, 46], as well as in the replica: in order to have a sensible 
Nf — ^ limit, the Goldstone manifold needs to be enlarged from SU(A^j) to U(iVj). It is then 
easy to see that the decoupling of the singlet field and the limit Nf ^ do not commute [37] . 
As a result, the chiral Lagrangian may contain all possible interactions involving the singlet, 
and to lowest order, the replica method has 



2Ne 2Nr 



where <I>o = (-F/2)Tr [— iln(C/)] and 17$ = exp(i6Ij^^jN^). Here, 1^^ is the identity matrix in 
the valence subspace and zero elsewhere. Obviously the couplings F and S need not be the 
same as in full QCD. In addition, new parameters related to axial singlet field, mo,a, have 
been introduced. In the SUSY formulation the first order chiral Lagrangian is the same, with 
the substitution Tr Str and U e U(iYi,|Ar„) (or, more precisely, U G Gi(Ar„|iV„) [47]). 

Since $0 is a singlet, there could in principle also be additional operators constructed with 
it in Eq. (5.1). They have, however, been shown to be suppressed by additional powers of 
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1/iVc [48], so we will neglect them in the following. 

Even though the low-energy Lagrangians of the replica and SUSY theories have quite 
different dynamical degrees of freedom, it is believed that in perturbation theory all com- 
putations concerning physical observables are equivalent [37]. In the e-rcgimc, however, a 
non-pcrturbative definition is needed for the zero-momentum integration. Traditionally this 
could only be achieved with the SUSY formulation, but there have been recent developments 
whereby it is argued that replica integrations can also be performed non-perturbatively, and 
agree with SUSY integrations [49] . Here we carry out the perturbative part of the computa- 
tion with the replica method, and return to the zero-mode integrations later on. 

Provided that only sectors of fixed topology are considered, the rules of the e-expansion 
are as in Eq. (3.6). The masslcss non-zero mode Goldstone propagator is also needed. In the 
replica case, it is given by Eq. (3.8), with 



However, as mentioned, our previous results were obtained with a completely general E(x), 
and therefore we know that Eqs. (4.1), (4.10) are independent of its form. 

5.2. Currents and weak operators 

Let us then consider the left-handed currents, in the replica formulation. Since the current 
follows from the Lagrangian, cf. Eq. (2.14), and the additional degree of freedom $0 is a 
flavour singlet, nothing changes with respect to the unquenched case at the present order: 



where the matrices are traceless and take non-zero values only in the valence sector. 

For the weak operators which do not directly follow from the Lagrangian, we have to be 
more careful. The general issue is whether there are more operators once the larger symmetry 
group of the quenched theory is considered. Clearly, for instance, one could attach the singlet 
field $0 to any operator. As mentioned, we assume such terms to be suppressed by 1/Nc 
and ignore them. Another trivial issue is that trace parts of {di^UU~^)ru vanish in the full 
theory but not in the quenched theory. However, there could in principle also be more drastic 
effects [24]. 

To be systematic about the operators appearing, let us recall the symmetries that are 
relevant. The weak operators we consider have indices corresponding to left-handed valence 
flavours only. Therefore they are singlets under the full right-handed symmetry group^, while 

®This excludes the problematic operators considered in ref. [24], ~ [Z7diag(/iv„, — /jv„) written here 

in the SUSY formulation. 




(5.2) 




(5.3) 



16 



they should have the correct symmetry properties under the left-handed valence subgroup. 
These requirements are sufficient to guarantee that the leading order 27-plet operator in the 

quenched theory is of the same type as in the unquenched case. 

Indeed, to get right-handed singlets under the full symmetry group, we are lead to the 
building blocks 

d^UU-\Ud^U~^ ~ 0(e2); UM, M^U'^ ~ 0{e^) , (5.4) 

which transform as fundamental (8> anti-fundamental under the valence subgroup SU{Nij)l, 
or SU(iV^|A^^)£,. We have also indicated the scalings of these operators in the e-rcgime. 
The operators can be trivially decomposed into a sum of 3 (g) 3* , 1 (g) 3* , 3 (g) 1 and 1 (g 1 
irreducible representations of the valence subgroup. To get a Lorentz invariant object with 
four flavour indices leads, as in the unquenched case, to a unique possibility up to and 
including 0{e^): 

[0^,]rsuv = ]F^(^^UU-') (di^UU-^) . (5.5) 

The reduction to irreducible representations follows from Appendix A. Only the 3 (g 3* 
components of d^UU~^ contribute to the 27-plet, because the 27-plet cannot appear in the 
tensor product of less than four fundamentals/anti- fundamentals. In other words, the oper- 
ator [OjuJ^su^ is symmetric in r s and u v and traceless in the SU(3)l subgroup, and 
zero if any of the indices lies outside of this subgroup. The weak Hamiltonian reads then 

= 2V2G^KdC{^52T"^'^'[0-]^„„. + •••} + H.C. , (5.6) 

where we have indicated that the quenched coupling, g^^^^'^^^'^ , does not need to be the same 
as the analogous one in the unquenched theory. 

In the case of the octets, on the other hand, the classification according to the valence 
group leads to ambiguities as discussed in [24]. 

Given that the currents have the same form as in the full theory, the quenched two-point 
function C"^''!"^"^^^'^ can now be obtained by setting Nf —>■ in the fixed-topology version 
of Eq. (4.1), as this result is independent of the part E[x) of the propagator. Alternatively, 
the result can be read from [44], by adding up the vector and axial- vector correlators and 
dividing by four: 

C*,u„o.=.(,„, ^ (-1VT"T')£ X {l + ^i- (ReTr. |Mf/„l)^^,_^ (5.7) 

where Tr^ denotes the trace over the valence subgroup, with = 3. Let us stress that the 
absence of E{x) guarantees that there is no dependence on the singlet couplings mg and a, 
which appear in the non-zero mode propagator of Eq. (5.2). Another interesting point to 
note is that the terms proportional to Nf of the unquenched result were largely responsible 
for the large corrections in realistic volumes, while they are absent now. This seems to imply 
that the e-expansion converges better in the quenched case. 
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Similarly, the result for the three-point function can be directly extracted from the fixed- 
topology version of Eq. (4.10): 



ol 



4L2 



F2 



Pi 



0«'00 



yV2 



[MUo, 



V 



(5.8) 



where we have omitted the indices from [A(^)]"^.j^„, [A(^)]"g„.y, given in Eqs. (4.8), (4.9). Again, 
the unphysical axial singlet couplings TOq, a do not appear at this order. 

An important point to stress is that the zero-mode integrals appearing in Eqs. (5.7), (5.8) 
are identical. Therefore, if we make the index choice of Eqs. (4.14), (4.15) and consider the 
ratio of Eq. (4.22), the results are the same in the unquenched and quenched cases: 



\ri 1 afe, quenched / \ 
2 QC'O-'^ , quenched (xg)C^^^ , quenched (yg) 



= 1 + 



F2 



/3i 



'0'>'00 



FV2 



V 



+ C»(e^) . 



(5.9) 



5.3. The quenched chiral condensate 

We end by discussing the actual value of the zero-mode integral in Eqs. (5.7), (5.8) for 
degenerate masses M = diag(m, m, m): 



^quenched — 



2^ .Tr.[^o + f/o- 



i/,Uo 



(5.10) 



This is just the quark condensate, and the vahic is well known, for A^.^ = 1 [43]. It is usually 
assumed that the outcome should not depend on > 1. What we wish to do here is to 
check this explicitly for A^^, = 2,3, by using the recent results from [49]: 

[J] . hm / det'^t/oe^^^^-I^^^"! - det [2^,(^ )] ^^.^^ ^ 

^/-o Juo 6 u(jv^) nf>i=i (m' - )«f ) 

where the NfXNf matrix Mj is Mj = diag(m-|- Ji, m-|- J^v^ , m, ...), and fii = {m+Ji)'EV, 
fi = mHV, together with 



i-i 



l(w), 



1,...,A^„ 

A-^ + l,...,2A'^ 



(5.12) 



The derivative of the logarithm of this partition function with respect to Ji, evaluated at 
Ji = 0, gives the required integral in Eq. (5.10) for any A^^,. The result indeed agrees for 
= 2,3 with the SUSY result obtained with U(l|l) [43]: 



^.quenched ^ Jj^^^^K^ifi) + I,+i{fi)K,_,{fi) 



+ - . 



(5.13) 
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For u ^ the leading behaviour of this function at small and large mass is the same as in 
the unquenched case at fixed topology, Eq. (4.21): 



6. Conclusions 

We have computed the two-point correlator of left-handed chiral charges, as well as the three- 
point correlator of two left-handed charges and one strangeness violating weak operator [13], 
in SU(3) chiral perturbation theory in a finite volume and close to the chiral limit, at ncxt- 
to-leading order in the e-expansion. The comparison of these observables with lattice data 
would in principle permit the extraction of the pion decay constant F, as well as the low- 
energy constant §27, involved in the A7 = 3/2 kaon decays and in the kaon mixing parameter 
Bk, with a minimal contamination from higher order corrections. Whether this will be 
numerically feasible is still an open question. 

We have also performed the same calculations in the quenched theory, using its replica 
formulation, and shown that these observables are only moderately affected. In particular, 
the ratio defined in Eq. (4.22) is not only independent of the quark masses sufficiently close 
to the chiral limit, but also unaffected by quenching at this order. Obviously the low-energy 
constants obtained with the quenched theory nevertheless differ from those in full QCD. 
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/X < 1 
Ai> 1 



(5.14) 
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Appendix A. SU(3) classification 

For completeness, we reiterate in this Appendix some essential aspects of the SU(3) classifi- 
cation of four quark operators. We follow the tensor method discussed, e.g., in [50]. 

The tensors we need to consider are of the form Orsuv, symmetric under {r s,u ^ v), 
and transforming under 3* (g) 3* (g) 3 (g) 3 of SU(3). We then define the projected operators 

^rsuv — (-^1 )rsuv;fsuv Ofsuvi (-'^•1) 
Orsuv — (-^2 ^rsuv;fsuv Ofgy^y, (■^•2) 

where cr = ±1. Here, with some redundancy in the symmetries of Pf , 

[Nf + 2a){Nf + a) 

+ + adsu^rf^vu^sv), (A. 4) 



Nf + 2(7 

where Nf = 3. 

It is easy to see that the antisymmetric tensor O^g^^ vanishes identically. The reason 
is that (as can be understood for instance by contracting with ^krs^iuv 

) it corresponds to a 

representation with dimension 8 just like RJ/^ defined in Eq. (A. 7), but all such representations 
have already been subtracted by the projection operator in Eq. (A. 4). 
Consequently, the reduction of a general operator Orsuv proceeds as 

1 



a=±l 



3(3 + a 



where O^guv transforms under the representation with dimension 27, and 



(A.5) 



^ Oli,i, (A.6) 

Here R^Js have the dimension 8, while S"^ arc singlets. 

Finally, let us note that in the chiral theory, i.e. if we replace Orsuv [Ow]rsuv^ 

so that 

[^w]rkuk ~ [^wV = 2^ [^tw]/feifc/' [''^u']™ ~ ~['^w]ru- (^--9) 
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